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We show that the transverse-momentum spectra of all hadrons measured at RHIC, 
including the hyperons, are described very well in a thermal model assuming the simul- 
taneous chemical and thermal freeze-outs. The model calculation takes into account all 
hadronic resonances and uses a simple parametrization of the freeze-out hypersurface. 

We present a simple model describing the p_i_-spectra of hadrons measured at RHIC at 
^/snn = 130 GeV. Our approach is a combination of the thermal model, used frequently in 
C\| ! the studies of the relative hadron yields with a model of the hydrodynamic expansion 
of matter at freeze-out. The main assumptions of the model [|10|-|T2|| are as follows: i) the 
chemical freeze-out and the thermal freeze-out occur simultaneously, which means that 
we neglect elastic rescattering after the chemical freeze-out, ii) all hadronic resonances 
are included in both the calculation of the hadron multiplicities and the spectra, and Hi) 
a simple form of the freeze-out hypersurface is proposed, which is a generalization of the 
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Bjorken model |J (see also EJ-QJ 



t 2 — x 2 — y 2 — z 2 = const. (1) 

The hydrodynamic flow on the freeze-out hypersurface (fl|) is taken in the form resembling 
the Hubble law, 

rf.= ^=*fi,£»,£V ( 2) 



r V 't'i't 

Recently, new arguments have been accumulated in favor of our first assumption. The 
measurements of the i^*(892) states by the STAR Collaboration j2D| indicate that either 
the daughter particles from the decay K*(892) — ► Kir do not rescatter or the expansion 
time between the chemical and thermal freeze-out is short (smaller than the K*(892) 
lifetime, r=4 fm/c). Moreover, the measured yield of if* (892) fits very well to the pattern 
obtained from the thermal analysis of the ratios of hadron abundances. This fact suggests 
again a short expansion time between the two freeze-outs. The assumption about the 
single freeze-out also solves the antibaryon puzzle pl |. Since the annihilation cross section 



for pp pairs is much larger than the elastic cross section, most of the protons would 
annihilate with antiprotons during the long way from the chemical to the thermal freeze- 
out. Such effect is not seen. In addition, let us mention that the single-freeze-out scenario 



* Supported in part by the Polish State Committee for Scientific Research, grant 2 P03B 09419. 



2 




0.5 1 1.5 2 2.5 3 3.5 



P ± [GeV] 

Figure 1. The p_i_-spectra at midrapidity of tt~, K~ ,p , and X*(892) in part (a), and of 
the hyperons A , H and f2 in part (b). The model calculation is compared to the PHENIX 
(filled symbols) and STAR (open symbols) most central data p0| , p3| -p0[ from Au + Au 
collisions at ^/snn = 130 GeV. Both the data and the theoretical curves are absolutely 
normalized (they include full feeding from the weak decays). 



is natural if the hadronization process occurs in such a way that neither elastic or inelastic 
processes are effective. An example here is the sudden- hadronization model of Ref. . 



Our model has two thermodynamic and two geometric (expansion) parameters. The 
two thermodynamic parameters, T = 165 MeV and \xb = 41 MeV, were obtained from 
the analysis of the ratios of the hadron multiplicities measured at RHIC ||. In this 
calculation the grand-canonical ensemble was used without the strangeness suppression 
factor (7 S =1). Since the particle ratios depend weakly on the centrality of the collision, 
we treat the thermodynamic parameters as the universal parameters (independent of 
centrality). The two geometric parameters are r of Eq. (|l|) and p max - The parameter 
Pmax determines the transverse size of the firecylinder at the freeze-out, 



p = yx 2 + y 2 < p max . (3) 
In the natural way, the values of r and p max depend on the considered centrality class of 



3 



events. For the minimum-bias data, which average over centralities, we find: r = 5.55 fm 
and p m ax — 4.50 fm, whereas for the most central collisions we find: r = 7.66 fm and 



Pmax = 6.69 fm HlQ| . The calculation of the spectra (and determination of the geometric 
parameters) is based on the standard Cooper- Frye formalism. The details of our method, 
especially of the technical problems concerning the treatment of the resonances, are given 
in the Appendix of Ref . [|TT|] . 

In Fig. 1 we show our results for the most central collisions. In the upper part (a) we 
show the spectra of pions, kaons, antiprotons, the (f) mesons, and the i\~*(892) mesons. In 
the lower part (b) we show the spectra of the hyperons A, H and Q. The model calculation 
agrees very well with the data. Note, e.g., the convex shape of the pion spectrum, crossing 
of the pion and the antiproton spectra at p± ~ 2 GeV, and the good reproducing of the 
fl spectrum. The good agreement between the model calculation and the data supports 
strongly the idea of thermalization of the hadronic matter produced at RHIC. Let us 



emphasize that the expansion parameters were fitted in Ref. |10| to the spectra of pions 



kaons and protons only. The spectra of other particles were calculated with the same 
values of the parameters, hence, they are predictions of our model. In view of this fact, 
the good agreement of the fl spectrum, predicted before the data were available, is highly 
non-trivial, especially in the context of the SPS results |I2]| . f\ 

A characteristic feature of our approach is a rather high decoupling temperature T ~ 
165 MeV. However, the (inverse) slope parameters corresponding to this temperature 
are lowered by the decays of the resonances [§]. This "cooling" of the spectrum by the 
decays of the resonances explains the difference between the high temperature of the 
chemical freeze-out and a smaller "apparent" temperature inferred from the shape of the 
spectra. We note that a similar high decoupling temperature has been found in the full 



hydrodynamic calculation of Ref. |31| , where also a complete set of hadronic resonances 
is employed. It remains a challenge to check whether our particular freeze-out conditions 
(shape of the freeze-out hypersurface and Hubble flow) may be obtained as the final stage 



of hydrodynamic evolution. First steps in this direction have been already made [p2 |. 

Let us make a few comments about the size of our geometric parameters. Translated to 
the measured HBT radii, R OVLt and -R S ide, they turn out to be too small. This problem can 
be circumvented by the inclusion of the excluded- volume corrections 0] which affect only 
the overall normalization of the spectra. If we rescale r and p max by about 30%, we obtain 
a satisfactory agreement with the HBT data. On the other hand, the ratio -R ou t/-R S ide 
is close to unity in our model, independently of the excluded-volume corrections. The 
approximate equality of these two radii follows in our model from the fact that the time 
extension of our system at freeze-out is much shorter than its space extension. 

In conclusion, we want to stress that a simple thermal model (with altogether four 
parameters) reproduces the transverse-momentum spectra of all hadrons which have been 
measured so far at RHIC. This fact brings strong evidence for thermalization of hadronic 
matter at RHIC and, possibly, indicates that a thermalized system of quarks and gluons 
was formed at the earlier stages of the collisions. 



We thank Boris Hippolyte for pointing out his talk J3U] with the experimental results 



2 The newest data on S production shown by Castillo during this conference have smaller overall normal- 
ization and also agree very well with our model calculation. 
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